Carrier transport in 2D graphene layers 
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Carrier transport in gated 2D graphene monolayers is theoretically considered in the presence of 
scattering by random charged impurity centers with density m. Excellent quantitative agreement 
is obtained (for carrier density n > 10 12 cm -2 ) with existing experimental data (Ref. Q], 0, @, EL H) ■ 
The conductivity scales linearly with n/ni in the theory, and shows extremely weak temperature 
dependence. The experimentally observed asymmetry between electron and hole conductivities is 
explained by the asymmetry in the charged impurity configuration in the presence of the gate voltage, 
while the high-density saturation of conductivity for the highest mobility samples is explained as a 
crossover between the long-range and the point scattering dominated regimes. We argue that the 
experimentally observed saturation of conductivity at low density arises from the charged impurity 



induced inhomogeneity in the graphene carrier density which becomes severe for n ; 



10 12 cm" 



PACS numbers: 81.05.Uw; 72.10.-d, 73.40.-c 

Recent experimental observation of the density de- 
pendent 2D carrier transport in single monolayers of 
graphene is an important, perhaps seminal, development 
in low dimensional electronic phenomena in nanostruc- 
tures [J, i, SHU. First, the gated 2D graphene 
systems could potentially become multifunctional high 
speed and high-power transistors, thus introducing a 
prospective paradigm shift in the micro- (and nano-) elec- 
tronics of the future. Second, carriers (both electrons 
and holes) have intriguing (and conceptually novel) lin- 
ear "Dirac-like" bare kinetic energy dispersion spectra in 
2D graphene monolayers due to the interesting honey- 
comb lattice structure of Carbon atoms (with 2 atoms 
on inequivalcnt sublattices per unit cell). Graphene is, 
in fact, a carbon nanotube rolled out into a 2D sheet, and 
the band structure induced carrier spectrum in graphene 
monolayers is thus linear and chiral with a four-fold 
ground state degeneracy arising from spin and valley, 
as well as a Berry phase term arising from the inher- 
ent sublattice symmetry which also restricts the carriers 
from back-scattering. This theoretical absence of back- 
scattering has led to the speculation that, as a matter of 
principle, carrier mobilities in 2D graphene monolayers 
could be extremely high even at room temperature, and 
indeed, already within the first six months since the be- 
ginning of the subject (i.e. six months since the fabrica- 
tion of the first gated 2D graphene samples with variable 
gate voltage tuned carrier density) mobilities as high as 
15000 cm 2 /Vs have been reported from IK to room tem- 
perature. These mobility values are comparable to (much 
higher than) the mobility in the best available Si MOS- 
FET samples at low (high) temperatures. It is, therefore, 
of considerable fundamental and technological interest to 
investigate the theoretical limit on the 2D graphene car- 
rier mobility. 



charged impurity centers, which is the most likely scatter- 
ing mechanism limiting graphene conductivity. We also 
calculate the effect of short-range scattering (e.g. by lat- 
tice defects), which may be present in some graphene 
samples. Quantitative agreement between our theory 
and existing graphene experimental transport data is 
a strong indication that the dominant carrier scatter- 
ing mechanism operational in 2D graphene monolayers 
is Coulomb scattering by random charged impurities lo- 
cated in the substrate near the interface between the 
graphene layer and the substrate. We estimate the typ- 
ical random charged impurity (effective) concentration 
to be approximately 10 12 cm~ 2 in currently available 
graphene samples, and suggest that reducing this impu- 
rity concentration to the 10 10 cm -2 density range - a dif- 
ficult but not an impossible materials fabrications task - 
should increase 2D graphene mobility to the extremely 
large value of ~ 1.5xl0 6 cm 2 /Vs even at high tempera- 
tures. We argue that the experimentally observed "low 
density" (< I0 12 cm -2 ) saturation of 2D graphene con- 
ductivity (i.e. the observed "minimum conductivity" at 
nominally zero gate voltage) also arises from the presence 
of these charged impurities which lead invariably to large- 
scale density inhomogeneities in the 2D carrier system 
with the 2D density fluctuations being larger than the 
average density at low carrier density (i.e. near zero gate 
voltage) where the system therefore breaks into a ran- 
dom network of 2D electrons and hole conducting pud- 
dles, producing a finite conductivity at zero gate voltage. 
Around zero gate voltage, the 2D graphene layer is thus 
not a simple (homogeneous) zero-gap semiconductor (as 
intrinsic graphene is expected to be), but a spatially in- 
homogeneous semi-metal with small random puddles of 
electron and hole 2D liquids depending on the details of 
the charged impurity configuration in the sample. 



In this Letter, we present a theory for 2D graphene car- 
rier transport taking into account scattering by random 



For our purpose, intrinsic graphene is essentially a 
zero-gap semiconductor (with the Fermi level Ep pre- 
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cisely at E = 0), with linear chiral carrier kinetic energy 
dispersion (arising from the known bulk band structure 
calculations) given by E — fvyk : where k is the 2D car- 
rier wavevector and j(= vp) is the constant (i.e. inde- 
pendent of carrier density) Fermi velocity. The intrin- 
sic situation (i.e. with zero gate voltage) has no free 
carriers at T — as a matter of principle (but, as a 
matter of practice the system breaks up into spatially 
inhomogeneous conducting puddles of 2D electron and 
hole droplets due to the potential fluctuations induced 
by the extrinsic random charged impurity centers which 
are invariably present in any real 2D graphene sample). 
The application of an external gate voltage (positive or 
negative) leads to free carriers (electrons or holes) in the 
system. We start by assuming the system to be a ho- 
mogeneous 2D carrier system of electrons (or holes) with 
a carrier density "n" induced by the external gate volt- 
age V g with n = K s V g /(4:irt), where k s is the dielectric 
constant of the substrate and t is the substrate thickness 
(i.e. the distance of the gate from the graphene layer). 
The chemical potential (or, equivalently the Fermi en- 
ergy Ep) at T = is given by Ep — hykp, where the 
2D Fermi wavevector kp depends on the carrier den- 
sity through kp — (Aim/ g^v) 1 / 2 with g s (g v ) being 
the spin (valley) degeneracy of graphene. (We choose 
Ss = 9v = 2 throughout this manuscript, although it is 
possible in some special situations for these degeneracies 
to be lifted.) 

For large carrier densities, where the system is homo- 
geneous, we expect the full Boltzmann transport theory 
developed below to be exact, hence providing an ac- 
curate determination of the charge impurity density rii 
by comparing our theory with experimental conductiv- 
ity data, since rii just sets the overall scale of the con- 
ductivity in the theory. This in turn provides an es- 
timate for the breakdown of our transport theory (typi- 
cally around n<rii) and the onset of the inhomogeneous 
puddles of electron and hole carriers. The experiments of 
Ref.0 show that the characteristic density for this break- 
down (i.e. the onset of inhomogeneity in our opinion) is 
n c ~ rii w 10 12 cm" 2 . An independent estimate (giving 
the same order of magnitude) can be made by considering 
the typical local density fluctuations induced by a single 
charge located inside the substrate at a distance d from 
the interface. This density fluctuation becomes global in 
the presence of multiple impurities which can be seen by 
considering multiple charged impurities distributed in a 
random Poisson manner in a 2D plane at a distance d 
above the graphene monolayer. One finds that the den- 
sity fluctuation Sn around the gate-induced average 2D 
carrier density n is given by Sn 2 = rii / (Sir d 2 ) confirming 
that transport close to the Dirac point is dominated by 
a spatially inhomogeneous carrier density. The transport 
through such 2D puddles is essentially a one dimensional 
random network for each of the four spin/ valley chan- 
nels providing a conductivity a ~ e /h, although a more 



complete theory is necessary to provide any quantitative 
comparison with experiments for n < n c ~ rii. 

We now proceed to describe in detail the microscopic 
transport properties at high carrier density using the 
Boltzmann transport theory Q. We calculate the mo- 
bility in the presence of randomly distributed Coulomb 
impurity charges near the surface with the electron- 
impurity interaction being screened by the 2D electron 
gas in the random phase approximation (RPA). The 
screened Coulomb scattering is the only important scat- 
tering mechanism in our calculation. There are addi- 
tional interface-scattering mechanisms unrelated to the 
Coulomb centers (e.g. surface roughness scattering), but 
such interface scattering is believed to be less important 
in the case of graphene (see below). 

We also neglect all phonon scattering effects mainly 
because both the reduced phase space imposed by chi- 
rality and the low temperatures cause phonon scattering 
to be negligible in the regime of current interest. Given 
that 2D graphene is essentially a weakly interacting sys- 
tem with effective r s ~ 0.75, a constant independent of 
carrier density (r s here is just the effective fine struc- 
ture constant), we expect our Boltzmann theory to be a 
quantitatively and qualitatively accurate description of 
graphene transport for all practical purposes. 

In Boltzmann theory the conductivity for graphene is 
given by a = (e 2 /h)(2Ep(r) /H) which comes from the 
massless chiral Dirac spectrum, where (r), the energy 
averaged finite temperature scattering time, is given by 



/& fc£fc r( efc )(-£) 



(1) 



where /(efc) is the Fermi distribution function, f(ek) — 
{l + cxp[(e^ — /i)]/fcsT}" 1 with /i(T, n) as the finite tem- 
perature chemical potential determined self-consistently. 

The energy dependent scattering time r(efe) for our 
model of randomly distributed impurity charge centers 
is given in the leading-order theory by 
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where is the concentration of the a-th kind of im- 
purity center, q = |k — k'|, 9 = # kk ' is the scatter- 
ing angle between the scattering in- and out- wave vec- 
tors k and k', e k = fi-7|k|, and v^(q,d) is the ma- 
trix elements of the scattering potential between an elec- 
tron and an impurity. For Coulomb interaction, we use 
y( c \q,d) = 2-7re 2 exp(— qd)/(kq) where d is the location 
of the charge impurity in the substrate measured from 
the interface and for short-range point defect scatterers, 
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v^(q,d) = vq, a constant. In Eq. [2l e(q) is the 2D finite 
temperature static RPA dielectric (screening) function 
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FIG. 1: Graphene conductivity limited by Coulomb scatter- 
ing calculated using different approximation schemes with di- 
electric constant k = 2.5. The complete screening approxi- 
mation (Ref. @) does not depend on dielectric constant and 
overestimates the conductivity, whereas using the unscreened 
dielectric function (also in Ref. @) would have conductivity 
less than 4e 2 /h for the entire range of gate voltages used in 
the experiment. RPA is the main approximation that we use 
in this work, and HA is the Hubbard approximation which 
incorporates local field corrections. 



where q s — Ae 2 kp / (hk-y) is the effective graphene 2D 
Thomas-Fermi wavevector and h — k(1 + irr s /2) is the 
effective dielectric constant. In calculating the conductiv- 
ity in Eq. [2j the factor ( 1 + cos &)/2 arises from the sublat- 
tice symmetry (related to Berry phase) which suppresses 
the backscattering contribution to resistivity. We can 
immediately observe from Eq. [21 that t -1 = r^ 1 + r" 1 , 
where the Coulomb scattering time r c ~ y/n domi- 
nates at low density over the short range scattering time 
t s ~ 1/s/ri. 

We emphasize that in order to get quantitative agree- 
ment with experiment, the full RPA dielectric function 
needs to be calculated. In Fig. [1] we show comparison 
between calculated conductivities in the RPA screening 
used in our calculation with other approximation schemes 
including the completely screened ( "strong screening" ) 



appropriate for graphene Q , given by 



if q < 2k F 




FIG. 2: The effect of remote scatters. Here d is the distance 
between the 2D graphene layer and the 2D impurity layer. 



approximation |8J valid only for r s ^> 1, which gives a 
larger value for the conductivity that is independent of 
the background lattice dielectric constant k. We note 
that for graphene r s ~ 0.75 making the strong screening 
approximation invalid. We also show results for Hub- 
bard approximation (HA) screening which includes local 
field corrections approximately. A recent theory Q has 
considered the temperature dependence of conductivity 
within the Thomas-Fermi approximation. 

Our main result is the quantitative agreement with ex- 
periments in the regime where the conductivity is linear 
in density. We note, however, that electron-hole asym- 
metry with a super-linear conductivity is ubiquitous in 
the experimental data. Changing the sign of bias volt- 
age has no effect on the concentration of charged impu- 
rities, but (and especially at high voltages) it changes 
the average distance between 2D graphene sheet and the 
impurity in the substrate, particularly since the random 
charged impurities could be both positive and negative 
in unequal numbers. In our model, we approximate that 
the impurities are confined in a 2D plane located at a 
distance d from the interface. We find that a small shift 
of d ~ 1 — 2 A is sufficient to explain the super-linear 
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FIG. 3: Graphene conductivity calculated using a combina- 
tion of short and long range scatterers. One finds that sub- 
linear conductivity at high density is likely to be seen in sam- 
ples with a small Coulomb impurity density and high mobility. 



behavior seen in the experiments. As shown in Fig. [21 
we notethat for fixed d, the difference between the actual 
conductivity and the linear value at d = is negligible 
at low density, n/rii ~ 1, but becomes important for 
n/rii > 3. 

Another feature of the experiments is that for very 
high mobility samples, one finds a sub-linear conductiv- 
ity [J]. Such high quality samples presumably have a 
small charge impurity concentration n% and it is therefore 
likely that point defects here play a more dominant role. 
Point defects (which would be dominant for either large 
carrier density n or low charge impurity concentration 
Hi) gives rise to a constant conductivity in contrast to 
charged impurity scattering which produces a conductiv- 
ity roughly linear in n/rii. Any realistic graphene sample 
is somewhere in the crossover between these two limits. 
Our formalism (see Eq.[2]) can include both effects, where 
zero range scatterers are treated with an effective point 
defect density of n p . 

We estimate that for most samples, we have n v jn% <SC 
1, while for the highest mobility samples, rip/fii ~ 0.2. 
Shown in Fig. [3] is the graphene conductivity calculated 
including both charge impurity (rii) and zero range point 
defect (rip) scattering for different ratios of the point scat- 
terer impurity density n p and the charge impurity den- 
sity Up. For small n p /rii we find the linear conductivity 
that is seen in most experiments. For large n p /n,i we 
see the flattening out of the conductivity curve (which in 
the literature [3J] been referred to as the sub-linear con- 
ductivity). We believe this high-density flattening of the 
graphene conductivity Q is a non-universal crossover be- 
havior arising from the competition between two kinds of 
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FIG. 4: Figure shows comparison of theory (Eq. [2} 
with experimental data of Ref. @ (up and down triangles, 
m = 2.3xl0 12 cm~ 2 ), Ref. (circles and squares, rn = 
3.4xl0 12 cm -2 ) and Ref. 3 (diamonds and crosses, rii = 
0.43xl0 12 cm~ 2 ) for electrons and holes respectively. Solid 
lines (from bottom to top) show the minimum conductivity 
of 4e 2 //i, theory for d=0, and d — 2 A. The inset shows Ref. @ 
data on a linear scale assuming that the impurity shifts by 
d — 2 A for positive voltage bias. 

scatterers, n p and rii. 

Finally, in Fig. [4] we show the comparison of our theory 
(with only charged impurity scattering) to recent exper- 
iments 0, 0, S| • We determine rn (which simply sets the 
overall scale) by fitting the data at high density to the 
result of Eq.[5]and observe very good agreement between 
all three experimental data sets and our transport theory. 
We note that the order of magnitude for m ~ 10 12 cm~ 2 
agrees with independent estimates (e.g. the low-density 
n < rii conductivity saturation around zero gate volt- 
age). The one-parameter (i.e. rij) quantitative agree- 
ment between our theory and three independent sets of 
experimental data is one of our most important findings. 

We emphasize that our Drude-Boltzmann semiclassical 
transport theory is valid only in the high density, n > n,, 
regime and therefore only the agreement between theory 
and experiment in the high density regime is meaningful. 
Indeed, our theory predicts that conductivity should van- 
ish linearly as carrier density goes to zero, whereas in re- 
ality, the graphene conductivity at low density (i.e. gate 
voltage), n < rii in our model, becomes approximately a 
density independent constant of the order e 2 /h. Our the- 
ory provides no quantitative or qualitative explanation 
for this low-density behavior, where we believe that the 
physics of strong density inhomogeneity (random pud- 
dles of 2D electron and holes) induced by the charged 
impurities in the insulating dielectric dominates, and our 
Drude-Boltzmann physics simply does not apply. We 
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cannot, however, rule out that the low-density regime is 
dominated by localization physics [l(| , or the Dirac cone 
physics of intrinsic graphene since our semiclassical 
theory does not apply in the n < rij low density regime. 
We note, however, that recent numerical work [a, Il2j] on 
quantum transport in the low density regime is consis- 
tent with the picture of a percolation transition between 
electron and hole carriers with the implication that the 
zero-carrier density intrinsic regime may simply be ex- 
perimentally inaccessible. However, what we have con- 
vincingly established in this work is that the high-density, 
n > rii, 2D graphene transport is dominated by impu- 
rity scattering which can be theoretically described by a 
microscopic Drude-Boltzmann model. 

In conclusion we have developed a detailed microscopic 
transport theory for 2D graphene layers. We have shown 
that charged impurities in the substrate are the dominant 
source of scattering. We have argued that the low den- 
sity regime is dominated by density fluctuations caused 
by these charge impurities and that a picture of inhomo- 
geneous puddles of conducting electrons and conducting 
holes is necessary to understand the finite value of con- 
ductivity near zero bias. Finally the Fermi temperature 
of 2D graphene being around 1300K for n ~ 10 12 cm" 2 , 
there is essentially no temperature dependence in 2D 
graphene conductivity in the 0—300K temperature range 
arising from charged impurity scattering, a fact that is in 
excellent agreement with experimental observation 13 1. 
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